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A new mechanism for spiral vortex nucleation in nongra- 
dient reaction diffusion systems is proposed. It involves two 
key ingredients: An Ising-Bloch type front bifurcation and 
an instability of a planar front to transverse perturbations. 
Vortex nucleation by this mechanism plays an important role 

^+in inducing a transition from labyrinthine patterns to spiral 

^^turbulence. 
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Labyrinthine patterns provide a beautiful example of 
pattern formation phenomena near equilibrium. Patterns 
of this form have been observed in various gradient sys- 
tems including garnet layers [1], ferrofluids [2] and block 
'copolymers [3]. Recently, Lee et al. found labyrinthine 



^-Npatterns in a bistable chemical reaction maintained far 

^—^from equilibrium [4]. Unlike equilibrium systems, non- 

^^gradient systems of this kind may undergo nonequilib- 

"^"rium Ising-Bloch transitions [5-7]. Such transitions or 

^^bifur cat ions render stationary (Ising) fronts unstable and 

^— Hgive rise to pairs of counter-propagating (Bloch) fronts. 

^r^[n this Letter we show the approach to an Ising-Bloch bi- 

I furcation may induce a transition from labyrinthine pat- 

H-^terns to "spiral turbulence", that is, a state of spatiotem- 

^poral disorder with repeated events of spiral vortex nu- 

'^leation and annihilation. 

Spiral turbulence has been found recently in reaction- 
diffusion models of cardiac tissues [8-10] and surface re- 
actions [11]. In both contexts wavetrain instabilities were 
found to precede the break up of spiral waves into new 
pairs of spiral waves. In this sense the resulting state of 
spiral turbulence is analogous to defect mediated turbu- 
lence in dissipative systems undergoing phase instabili- 
ties [12]. Here we report a different mechanism for spiral 
vortex nucleation that involves two generic front insta- 
bilities: An Ising-Bloch front bifurcation, and an insta- 
bility to transverse perturbations [13,14]. The front bi- 
furcation provides the system with two coexisting Bloch 
fronts. The transverse instability provides the driving 
force needed to convert a segment of a Bloch front prop- 
agating in one direction into a Bloch front propagating in 
the opposite direction. This process nucleates a pair of 
spiral vortices, one at each end of the segment. Such nu- 
cleation events can be realized by approaching the front 



bifurcation point from the regime of labyrinthine pat- 
terns, for the latter already imply the existence of a trans- 
verse instability. 

We study spiral vortex nucleation using a nongradient 
reaction-diffusion model of the FitzHugh-Nagumo type: 
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where u and v are two scalar real fields. The model con- 
tains four parameters: e > 0, the ratio between the time 
scales associated with the two fields; ^ > 0, the ratio 
between the two diffusion constants; and the two param- 
eters, ai > and qq. We choose qq and ai so that the sys- 
tem (1) has two stable uniform steady states: A "down" 
state (i/_,i;_) and an "up" state (i/+,i;+). Note that 
for ao = the up and down states have the symmetry 
(i/_,i;_) = — (i/+,i;+). This model has the interesting 
property that for e ^ 1 it can be reduced to a nonlocal 
gradient system by eliminating v adiabatically [15]. Sim- 
ilar nonlocal models have been studied in the context 
of phase separating systems [16,17]. The regime e ^ 1 
has been studied most recently by Petrich and Goldstein 
who derived heuristically a nonlocal equation of motion 
for the front interface separating the up and down states 

[IS]- 

The Ising-Bloch front bifurcation for the system (1) 
has been studied in Refs. [7,19]. In the regime 6/e ^ 1 
and for ai sufficiently large, the speed c of a planar front 
connecting the up state at a? = — cxd to the down state at 
X = DO satisfies 
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In the 



where r]^ = eS, q^ = ai -\- 1/2 and Coo = -j=^ 
symmetric case (ao = 0) a graph of c = c(r]) yields a 
pitchfork bifurcation diagram. For ri> ttjc — ^ 3 , c = 
is the only solution, representing a stationary Ising front. 
For rj < ?7c, two additional solution branches appear, c = 
:^2qy^r]'^ — rf)^ representing counter propagating Bloch 
fronts. Structurally, a Bloch front differs from an Ising 
front in that the v field is displaced with respect to the 
u field. The magnitude and the sign of the displacement 



determine the speed and the direction of propagation, 
respectively. A useful rule to remember is that the v 
field always lags behind the u field. 

A typical form of the bifurcation diagram for the non- 
symmetric case is shown in Fig. 1. Front multiplicity 
arises through a saddle-node bifurcation at a critical r] 
value, r]rXao)- We will refer to the single front solution 
that exists for 77 > 77^ as an Ising front and to the two sta- 
ble front solutions for 77 < 77^ as Bloch fronts. With this 
convention the Ising front and one of the Bloch fronts lie 
on the same solution branch. Note that the Ising front is 
not stationary. 
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FIG. 1. A front bifurcation diagram for the nonsym- 
metric case, obtained from equation (2). Parameters used: 
ao = —0.1, ai = 2. 

To study the transverse instabilities of the various front 
solutions we consider a weakly curved front. Introducing 
an orthogonal coordinate system (s, r) moving with the 
front, the equations (1) assume the leading order form 

,3 
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The coordinate s parametrizes the front and Cr{s^t) and 
K(s,t) are, respectively, its normal speed and curvature. 
Both Cr and k may vary weakly along the front and evolve 
slowly in time. Multiplying (36) by the factor A(s,t) = 
(cr + K)/(cr + ^k) the resulting system can be interpreted 
as equations for a planar front propagating at speed c^-hK 
along the r axis in an efi'ective medium characterized by 
the parameters e = eA and 6 = ^A [20]. Using equation 
(2) with c replaced hj c^ -\- k and 77 by 77 = 77A we obtain 
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Equation (4) can be solved for c^ in terms of k. Since k 
is small we expand c^ = cq + c?k -h C^(k^), where 00(77) is 



the speed of a planar front satisfying (2), and evaluate 
the prefactor d. This yields 
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Setting d = for a given branch of cq gives the parameter 
values at which the corresponding planar front becomes 
unstable to transverse modulations. For the symmetric 
case (ao = 0) the Ising and Bloch fronts become unsta- 
ble to transverse modulations when S > Sj(e) = g^^c 
and S > Sb(c) = ^-^3 .- , respectively. The transverse 

instability lines for the nonsymmetric case are shown in 
Fig. 2. Notice that the lines corresponding to the two 
Bloch fronts, denoted by S^, are not degenerate as in the 
symmetric case. Also shown in this figure is the front 
bifurcation line, 77 = 77^ which we write as ^ = ^i^(6). For 
the symmetric case Sp^e) = gT- 



^ 



6 
4 

2 


'5 = <5B*{e)^'^.,^^_^ 


— 1/ 




6 = 6.{e) 








0.0 



0.1 



FIG. 2. The front bifurcation line (thick) and the trans- 
verse instability lines (solid thin) for the nonsymmetric case, 
obtained from equations (2) and (5) respectively. The dashed 
line represents the transverse instability of the unstable front 
branch (dashed line in Fig. 1). Parameters used: ao = —0.1, 
ai = 2. 

The front bifurcation line and the transverse instabil- 
ity lines divide the e — 6 parameter plane into regions 
where difi'erent pattern behaviors are expected. In the 
Ising regime, 6 > ^i^(6), and below the transverse insta- 
bility line of Ising fronts, 6 < ^/(e), we find stationary 
stripe patterns or, for 6 suflEiciently small, no patterns at 
all. In the Bloch regime, and below the transverse insta- 



bility lines, S 



^S' 



traveling waves, including rotating 



spiral waves, prevail [7]. We emphasize though that the 
bifurcation and the transverse instability lines were de- 
rived for planar fronts and that curvature may shift these 
lines. In particular, it is possible to have two appropri- 
ately curved fronts propagating in opposite directions on 
the Ising side of the front bifurcation line. 



Fig. 3 shows the time evolution of a stripe of up state 
domain well inside the Ising regime and above the trans- 
verse instability line, obtained by numerical integration 
of (1). The stripe is transversally perturbed in its mid- 
dle part and at first a meandering stripe forms. At later 
stages, fingering and tip splitting processes take place 
until a stationary labyrinthine pattern fills up the whole 
system. An important characteristic of the final pattern 
is that the up-state domain is connected; domain fission 
is avoided because of repulsive front interactions. This 
behavior is very similar to that observed by Lee et aL in 
the bistable chemical recation [4] and to the simulations 
of Petrich and Goldstein [18] on their nonlocal interface 
model. We note that the value of e used in producing 
Fig. 3 was quite small (e = 0.05). This suggests the 
gradient nature of (1) for e ^ 1 persists for smaller e val- 
ues as well, extending the validity range of the nonlocal 
interface model. 
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FIG. 3. The development of a labyrinthine pattern from a 
single stripe in the Ising front regime. The light and dark 
regions correspond to the downand up states, respectively. 
The frames a,b,c,d pertain to times t = 100,950,1900,5000. 
Parameters used: ao = —0.1, ai = 2, e = 0.05, ^ = 4. 

SuflEiciently close to the front bifurcation line, however, 
nongradient efi'ects become important. Fig. 4 shows the 
time evolution of a stripe in the vicinity of that line. With 
this parameter choice both the leading and the trailing 
fronts of the stripe (Bloch fronts of difi'erent kinds) are 
unstable to transverse perturbations, but the instability 
of the leading front is stronger. In contrast to the time 
evolution shown in Fig. 3 the stripe breaks into disjoint 
pieces which subsequently develop into a complex spa- 
tiotemporal pattern. We have continued the simulation 
up to ^ = 5000 without finding any qualitative change. 



The number of vortices appears to fiuctuate between 30 
to 70. We point out, however, that equations (1) always 
have the two simple attractors, (i/_,i;_) and (i/+,i;+), 
and, in principle, further evolution might culminate in 
one of them. 




FIG. 4. The development of spiral turbulence from a single 
stripe close to the front bifurcation line. The light and dark 
regions correspond, respectively, to the down and up states. 
The frames a,b,c,d pertain to times i = 30,210,1180,1800. 
Parameters used: ao = —0.1, ai = 2,e = 0.014, 8 = 2.8. 

Fig. 5 provides a closer look at the processes involved 
in Fig. 4. Shown in this figure are the contours i/ = 
(thick line) and v — {) (thin line) for four consecutive time 
frames. The frames show a dent in a Bloch front prop- 
agating to the right that turns into a Bloch front prop- 
agating to the left. This process is accompanied by the 
formation of a vortex pair appearing as crossing points of 
the two zero contours. The growing dent eventually cuts 
the up-state domain into two disjoint pieces in a process 
involving down-state domain fission and front reconnec- 
tions. Domain fission is made possible at this parameter 
regime because of higher front speeds. 

The rate of vortex nucleation depends not only on the 
strength of the transverse instability but also on the oc- 
casional formation of highly curved front regions. Such 
regions are most often created by the fission of two up- 
state domains, yielding cusp like structures as shown in 
Fig. 5a. This geometry facilitates the difi'usive accumu- 
lation of i; at a small region, thereby inducing a local 
front transition and consequently vortex nucleation. 



The simulations shown in Figs. 3 to 5 were obtained by 
integrating the system (1) using an implicit finite differ- 
ence scheme for the labyrinthine patterns and an explicit 
scheme for the turbulent patterns. A 400 x 400 grid has 
been used with one grid point per unit length. This re- 
sulted in about 6 grid points across the front. Doubling 
the resolution slightly changed the location of the front 
bifurcation line but yielded the same qualitative results. 
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FIG. 5. Vortex nucleation, developing from a cusp like 
front shape, followed by domain fission and front reconnec- 
tions. The thick (thin) line is a m = (t? = 0) contour. The 
light regions correspond to the down state and the shaded 
regions to the up state. The frames a,b,c,d pertain to times 
t = 310, 390, 460, 500. Parameters are as in Fig. 4. 

We have shown here that the coupling between an 
Ising-Bloch bifurcation and a transverse front instability 
in a doubly-diffusive FitzHugh-Nagumo model may lead 
to spontaneous nucleation of spiral vortices. These nucle- 
ation events, together with the capability of domains to 
merge, can induce a transition to spiral turbulence. Sim- 
ilar behavior is expected to be found in other reaction 
diffusion systems exhibiting these two front instabilities. 
A possible candidate is the bistable reaction studied by 
Lee et aL [4]. 
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